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The compaction dynamics of a granular media subject to a sequence of vertical taps made of fluid 
pulses is investigated via Molecular Dynamics simulations. Our study focuses on three different 
levels: macroscopic (volume fraction), mesoscopic (Voronoi' volumes, force distributions) and micro- 
scopic (grain displacements) . We show that the compaction process has many characteristics which 
are reminiscent of the slow dynamics of glass forming systems, as previously suggested. For instance 
the mean volume fraction slowly increases in time and approaches a stationary value following a 
stretched exponential law, and the associated compaction time diverges as the tapping intensity 
decreases. The study of microscopic quantities also put in evidence the existence of analogies with 
the dynamics of glass formers, as the existence of dynamical heterogeneities and spatially correlated 
motion of grains; however it also shows that there are important qualitative differences, as for in- 
stance in the role of the cage effect. Correlations between geometry and dynamics of the system at 
the grain level are put in evidence by comparing a particle Voronoi' volume and its displacement in 
a single tap. 



I. INTRODUCTION 



When subject to vertical vibrations granular materials 
can produce a variety of distinct phenomena, depend- 
ing both on the driving parameters and on the container 
properties. A great deal of interest has been recently 
raised by the process of compaction under successive ver- 
tical taps [l|,EBi,ll- Experiments [H, i, Hi show that 
the density of a column of grains submitted to vertical 
tapping increases slowly in time with a law well fitted by a 
stretched exponential [l| or logarithmic function Q , and 
that the characteristic compaction time grows abruptly 
as the driving intensity decreases. These observations 
have suggested an analogy with glass-forming systems, 
where the relaxation time diverges as the temperature is 
decreased; the analogy is corroborated by the fact that 
concepts like frustration and free volume, which are com- 
monly used to explain the slow dynamics of supercooled 
liquids and other thermal systems, do also provide an 
insight in the physics of powder compaction [1, @, 0] ■ 

However, contrary to supercooled liquids, granular ma- 
terials are non-thermal systems as their typical energy 
scale, the energy required to rise a grain of its own diam- 
eter against gravity, is orders of magnitude larger than 
the thermal energy /cbT: each granular pack is in a me- 
chanically stable state, which last as long as there is no 
external perturbation. Therefore even though the dy- 
namics induced by a sequence of vertical taps becomes 
slower and slower as the tapping intensity decreases, im- 
portant qualitative differences may exist between slow 
granular dynamics and glassy dynamics. 
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Here we investigate analogies and differences between 
granular dynamics and glassy dynamics by performing 
molecular dynamics (MD) simulations of a granular sys- 
tem subject to a sequence of vertical taps where, in or- 
der to explore a wide range of volume fractions, the 
system is tapped via flow pulses as in the experiment 
of Schroter et. al. [5]. We investigate the evolution 
of macroscopic quantities (volume fraction), mesoscopic 
quantities (Voronoi volumes, force distributions) and mi- 
croscopic quantities (grain displacements), and we dis- 
cuss how various static properties of a granular pack 
change during compaction We found severals analo- 
gies between the compation of granular media, and the 
slow dynamics of glass forming systems [1, [13, [Hi, as 
the divergence of the relaxation times, dynamical het- 
erogeneities and spatially correlated motion; but we also 
show evidences of important qualitative differences, as in 
the role of the cage motion, which are put in evidence by 
the study of particle trajectories. 

The paper is organized as follows. Sec. [Til presents the 
numerical model used. Then compation dynamics, inves- 
tigated via the study of the time dependence of the vol- 
ume fraction, and of the diffusion coefficient at station- 
arity, is discussed in Sec. lHIl The evolution of structural 
properties of a compacting granular pack, the radial dis- 
tribution function, the distribution of Voronoi volumes, 
and the distribution of interparticle forces, is presented 
in Sec. [TV] Sec. |V| discusses the compaction dynamics 
at a grain level, showing the existence of dynamical het- 
erogeneities and of spatially correlated motion of grains, 
and putting in evidence qualitative differences in the par- 
ticle trajectories of compacting granular media and glass 
formers. Sec. |V| also presents a connection between geo- 
metrical (Voronoi volumes) and dynamical (particle dis- 
placements) properties of the system. Finally, a conclu- 
sion summarizes the main results and perspectives. 
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II. NUMERICAL MODEL 

We run Molecular Dynamics (MD) simulations of iV = 
1600 monodisperse spherical grains of diameter d — 1cm 
and mass m = Ig. Grains, under gravity, are confined in 
a box with a square basis of length L — 10cm, with pe- 
riodic boundary conditions in the horizontal directions. 
The bottom of the box is made of other immobile, ran- 
domly displaced, grains (to prevent crystallization) . Sim- 
ulations with 4 time more particles in a system with a 
square basis of length 2L give the same results, so we 
exclude the presence of finite size effect in our system. 

Two grains in contact interact via a normal and a tan- 
gential force. The former is given by the spring-dashpot 
model, while the latter is implemented by keeping track 
of the elastic shear displacement throughout the lifetime 
of a contact [l^]- The model is the one described in [l^ 
with a restitution coefficient, e = 0.8. We use the linear 
model insted of the more realistic Hertzian model [l3| 
as this latter is characterized by a coefficient of restitu- 
tion which goes to zero as the relative velocity of the 
impacting particles decreases [l^ . This feature makes 
computationally expensive the simulation of a granular 
system reaching a mechanically stable state. 

The system is immersed in a fluid and, starting from a 
random configuration, it is subject to a dynamics made 
of a sequence of flow pulses where the fluid flows through 
the grains (see FiglT]), as in the experiment of Ref.jH]- In 
a single pulse the flow velocity, directed against gravity, 
is y > for a time tq; then the fluid comes to rest. We 
model the fluid-grain interaction as in Ref.s [l^ [l^ via 
a viscous force proportional to the fluid grain relative 
velocity: Ffg — —A{v — V) where v is the grain and 
V = (0, 0, V) is the fluid velocity. The prefactor A = 
7(1 — $;)~3-65 jg dependent on the local packing fraction, 
<i>; , in a cube of side length 3d around the grain, and the 
constant is 7 = 1 Ns/cm [l6j . 

During each pulse, grains are fluidized and then come 
to rest under the effect of gravity g. The system is con- 
sidered to be still when the kinetic energy per grains is 
below 10~^mgd. All measures below are recorded when 
the pack is at rest. 

The dynamics of dry granular media subject to vertical 
vertical vibrations is determined by two parameters, the 
amplitude and the frequency of vibrations. In the system 
we are investigating here there are also two parameters, 
the tap duration tq and the fluid velocity V . 



III. DYNAMICS 

A. Volume fraction 

When subject to a sequence of flow pulses a granu- 
lar system compactifies (or expands) until it reaches a 
stationary state which depends on the driving parame- 
ters [11. Fig. [2] shows the evolution of the volume fraction 
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FIG. 1: (Color online) We study the compaction process of 
a granular media subject to flow pulses. During each pulse 
(tap), of duration ro, fluid flows with velocity V trough the 
granular media. Before applying a flow pulse we wait until 
the granular media comes to a rest. 
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FIG. 2: (Color online) Temporal evolution of the mean vol- 
ume fraction "I> for ro — 0.03 s and the reported values of 
fluid velocities V. The data obtained via simulations of a 
larger systems (4 times more particles) evidence the absence 
of finite size effects. Dashed lines are fit to a stretched expo- 
nential law. 



$ (measured in the bulk of the system) with the num- 
ber of flow pulses for systems subject to a tap dynamics 
with tap length tq = 0.03s and various values of the fluid 
velocity V. Each curve is obtained by averaging over 
32 independent realizations. Similar curves are obtained 
with different values of tq. 

The time evolution of the volume fraction is well de- 
scribed by a stretched exponential law. 



Ht) = $00 - (*oo - $0) exp [-{t/ry 



(1) 



with c ~ 1. This is in agreement with experimental re- 
sults of Philippe et al. [ll.l9l.IT7t. which have investigated 
the relaxation dynamics of dry granular media subject to 
vertical taps in a system with a height to width ratio sim- 
ilar to ours. On the contrary Nowak et. al. [3, 4] have 
investigated a system with a larger height to width ratio, 
finding logarithmic compaction. 

As in previous experiments of both vibrated [ll.l3l.[6l[T8j 
and fluidized [5| granular systems, when the tapping in- 
tensity decreases the system compactifies more. This is 
clearly illustrated in Fig. [31 where we show the depen- 
dence of the volume fraction reached at stationarity on 
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FIG. 3: (Color online) Dependence of the volume fraction 
reached at stationarity by a system subject to a sequence of 
flow pulses of length to on the fluid velocity V. Lines are a 
guide to the eye. 

the fluid velocity V for various values of the tap length 
To . As one could have expected the final stationary state 
depends both on V and tq. However it is possible to 
show numerically [l^ that the final state can be charac- 
terized by one thermodynamical parameter, supporting 
the idea of a statistical mechanics description of granular 
media at rest 0, H [m nil originally proposed by S.F. 
Edwards 

B. Compaction time 

The relaxation time t which appears in Eq.[l]is a mea- 
sure of the number of fiow pulses required by the system 
to reach stationarity. The experiments of Ref. Q have 
investigated a range of parameters V, tq where the system 
reaches a stationary state after few flow pulses. On the 
contrary, the experiments of Bideau et. al. [l| have inves- 
tigated a range of parameters where the compaction dy- 
namics of the system is glassy like. They report a relax- 
ation time following an Arrhenius behavior with the in- 
verse maximum acceleration of the pack, r oc exp (r~^). 
In the system under investigation here the relaxation 
time also evidences the existence of a glassy dynamics 
of the system, as it diverges with a power law with de- 
creasing fluid velocity, 

T OC (2) 

with P ~ 1.17 ± 0.04, as shown in Fig. |4] (for data with 
tap length tq ~ 0.03 s). The same behavior is observed 
for different values of the tap length tq. 

C. Stationary dynamics 

After having applied a long sequence of pulses up to 
reach the stationary state, we have computed the mean 



FIG. 4: The relaxation time r increases with a power law 
(Eq. [2]) as the fluid velocity decreases. 

square displacement, 

1 ^ 

i 

with ri{t) position of grain i after t taps, and tw waiting 
time which depends on the driving conditions {t^ — 3t). 
From the mean square displacement, shown in Fig. [51 
we have extracted the diffusion coefficient D ((r^(t)) (x 
Dt) which decreases with the fluid velocity as a power 
law, D cx with the same exponent observed for the 
dependence of the relaxation time on the fluid velocity 
(Fig. [5]). This suggest the existence of relation r cx Z?^^ 
relating the compaction time and the diffusion coefficient 
at stationarity, which is illustrated in the inset of Fig. [5| 
(lower panel). 

The mean square displacement evidences the absence 
of a subdiffusive regime in the slow dynamics of granular 
media subject to flow pulses, a signature of the cage ef- 
fect in supercooled liquids. In Sec. |V] we will show that 
particles may be constrained in a cage, but that the es- 
caping time is too small (few taps) in order to affect the 
mean square displacement. 

IV. STRUCTURE EVOLUTION 

A. Radial distribution function 

The radial distribution function g{r) is the probability 
distribution of flnding the center of a particle in a given 
position at a distance r from a reference sphere. Since it 
contains information about long range interparticle cor- 
relations, it is a common tool in the characterization of 
packing structures. Here we study how the radial dis- 
tribution functions, which we normalize as usual is such 
a way that g{r) — > 1 for r — s- cx), evolves during com- 
paction. 
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FIG. 5: (Color online) Left panel: mean square displacement 
(in the stationary state) for to = 0.03 s and fluid velocities 
(from top to bottom) V = 2.00, 1.25,0.8,0.40,0.15 and 0.10 
cm/s. Right panel: the inverse diffusion coefficient (measured 
for those points in which the diffusion regime is attained) di- 
verges as a power law, — aV^^ as the fluid velocity V 
decreases. Inset: Stokes-Einstein relation between the com- 
paction time and the diffusion coefficient at stationarity. 



Fig. [6] (main frame) shows g(r) at different times dur- 
ing compaction (tq = 0.03 s, V = 0.2 cm/s); similar 
results are found with different values of tq and V. The 
first strong peak at r = d = 1 cm corresponds to the 
high probability of a having a neighbor in contact. This 
peak characterizes all dense systems of hard particles, as 
a consequence of their impenetrability. In a granular me- 
dia at rest under gravity it is enhanced by the fact that, 
in order for the system to be stable, each grain must con- 
tact other grains. The following two maxima, as shown 
in the inset, appears at r = \/3d and r = 2c?. Both 
of them increase with the volume fraction indicating an 
increasing organization of the packing (not necessarily 
related to the formation of ordered structures [1^). A 
similar dependence of the secondary peaks of the radial 
distribution function on the volume fractijan has been ob- 
served experimentally by T. Aste at al. |25| , which have 
investigated via x-ray tomography packs with different 
densities, and numerically by L.E. Silbert at al. p^ . 
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FIG. 6: (Color online) Radial distribution function for pack- 
ing of volume fraction <I> — 0.576, 0.596 and 0.616. As the 
volume fraction increases we observe a small increment of the 
peaks at \/3d and 2d, with d — 1cm grain diameter. 



B. Voronoi tessellation 

In a granular pack of monodispere spheres of volume 
fraction $ the mean volume occupied by a particle is 
Vp/^, where Vp is the volume of a particle. When the 
system is in a disordered state there will be both parti- 
cles occupying a larger volume, and particles occupying a 
smaller one. It is therefore instructive to investigate what 
is the probability that a given particle occupies a volume 
V, and how this probability changes during compaction. 

To this end one has to operatively define what is the 
volume occupied by a particle: by using the Voronoi tes- 
sellation (as in [25|) we define the volume Vi occupied 
by particle i as the volume of the convex polyhedron 
which contains all points closer to particle i than to any 
other particle. Fig.jTja shows the distribution P{v) of the 
Voronoi volumes of a system tapped with tq — 0.03 s and 
V = 0.2 cm/s, after 1,10,100 and 300 taps. These are 
slightly asymmetric distributions with exponential tails 
(Fig. |7|a, inset). The asymmetry is a standard feature 
of the Voronoi distribution: for an ideal gas in one di- 
mension, for instance, P{v) cx vp^(l — pY , where the 
Voronoi volume w of a given atom is half of the distance 
between its left and right nearest neighbors. As the sys- 
tem compactifies both the mean value (w) and the stan- 
dard deviation of the distribution decrease. However 
the distribution retains its functional form: when ayP{v) 
is plotted versus {v — (v))/ all of the different curves 
scale on the same master curve, as shown in Fig. |7|b. 
This scaling suggests the existence of a single geometri- 
cal structure of the system, only specified by its volume 
fraction, and tell us that there are no dramatic structural 
changes during compaction. 

The same scaling of the distribution of Voronoi vol- 
umes has been observed by F.W. Starr tt al. [26| in MD 
simulations of a glass-forming polymer melt, and is ver- 
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FIG. 7: (Color online) (a) Distribution of the Voronoi volumes 
in the sample after 1, 10, 100 and 300 taps. As the system 
compactifies the standard deviation and the mean value {v) 
of the distribution decrease, (b) Scaling of the distributions 
of the Voronoi volumes shown in Panel (a). Same symbols are 
used, (c) The same scaling has been found in MD simulations 
of a model of glass-formes by Starr et al. JM, and is verified 
by the experimental data of Aste et al. [2g|. However data 
from different sources do not scale on the same master curve. 



ified by the experimental data on granular packs of T. 
Aste et al. [1^. However the data from these different 
sources do not scale on the same muster curve, as shown 
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FIG. 8: (Color online) The distribution of interparticle normal 
forces decay exponentially at high forces, it appears to be 
insensitive on the packing fraction of the system. 

in Fig. [71c. The discrepancy which may be due to the 
fact that Ref. |26;] investigates a thermal system, while 
Ref. investigates a polydisperse system, needs further 
investigation. 



C. Force distribution 

Since the work of Mueth et. al. j27^ interparticle force 
distributions have become a standard tool for the charac- 
terization of granular packs. There are now experiments, 
numerical simulations and theories (see [1^ and refer- 
ences therein) finding an exponential decay at high forces. 
This exponential decay is the signal of an heterogeneous 
structure of the system: while most of the interparticle 
normal forces have magnitude close to the mean value, 
there are also interparticle normal forces of much higher 
magnitude. Moreover, these high forces have been shown 
to be spatially related, giving rise to the well know force 
chains. 

We have studied the evolution of the probability distri- 
bution of normal forces during compaction. In order to 
avoid the effect of gravity (due to the absence of vertical 
walls the mean vertical stress depends on the depth) we 
have computed the probability distribution of the nor- 
mal forces between grains contacting in a point at a 
height z enclosed in a thin horizontal slice (z > 8 cm 
and z < 10 cm). Similar results are observed when se- 
lecting different horizontal slices of our system. Fig. [5| 
shows the force probability distribution after 1, 10, 100 
and 300 taps (corresponding to volume fractions in the 
range 0.576 — 0.616). Since the distributions collapse on 
the same curve normal forces appears not to be affected 
by the density of the system (in the range we have in- 
vestigated). We conclude that the force distribution is 
rather insensitive to the density of the granular pack, as 
also observed in Ref. [s^ and in [s^j (where it is show 
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FIG. 9; (Color online) The compressional force C of a particle 
decreases with its Vornoi' volume. 

that forces do not couple with the density for thermal 
systems near the glass transition as well). It should be 
noted, however, that the use of the more realistic Herzian 
(instead of the linear) grain-grain interaction model (see 
Sec. ini may change the properties of the force distribu- 
tion. 



FIG. 10: (Color online) Typical trajectory of a particle during 
compation, followed for 300 taps during compaction (V — 0.2 
cm/s. To = 0.03 s). Circles illustrate that sometimes a particle 
is confined in cages formed by its neighbors. 



reasons. However a reliable estimate of a is difficult to 
obtain as Voronoi volumes vary in a small range. 



V. GRAIN MOTION 



D. Force- Volume correlations 

The relation between geometrical structure of the 
packing and interparticle forces has been investigated in a 
number of previous works (see [11] and reference therein). 
Here we report on a relation between the Voronoi volume 
associated to a particle and the forces acting on it. To 
this end we define the compressional force Ct acting on 
particle as: 

a = (4) 

where \fij\ is the normal force of interaction between 
particles i and j. Ci measures how much particle i is 
compressed as the pressure acting on it is Ci/Si, where 
Si = 47r(d/2)^ is the particle surface. Fig. [3] shows that 
the compressional force decreases with the Voronoi vol- 
ume. 

A simple explanation of the decreasing of the compres- 
sional force with the Voronoi volume can be obtained via 
the following argument. Consider two contacting parti- 
cles at a distance I. Their interparticle force decreases 
with I and particularly in our case, due to the computa- 
tional model used, / = k{d — I), where d is the diameter 
of a particle. On the other hand the Voronoi volume v of 
one of these particles increases with I. Assuming v oc I", 
one obtains / oc {d—v^^" ), i.e. a decrease of the compres- 
sional force with the Voronoi volume. The data of Fig. [5] 
are consistent with a = 3, as expected for dimensional 



In this section we investigate how a granular pack 
moves during a single tap. We consider a system subject 
to a tap dynamics with tq = 0.03 s and V = 0.2 cm/s, 
and compare the state reached after the application of n 
taps, with the state reached after the application of one 
more tap (with n=l, 10, 100, 300). 

Particularly we investigate the distribution of particle 
displacements, the heterogeneity of particle motion, and 
the correlation between a particle displacement and both 
its Voronoi volume and the compressive force acting on 
it. 



A. Particle trajectory and cage motion 

Fig. rroi shows a typical trajectory of a particle during 
compaction. As observed in colloidal systems the trajec- 
tory is characterized by period of time in which the par- 
ticle is confined in cages formed by its neighbors. Cage 
motion has also been observed before in experiment of 
granular materials subject to continuous vibrations (29j 
and to shear ^30, 31, 32]. 

The typical linear size of the cage is roughly 0.05d, 
where d = 1 cm is the diameter of the particles. A sim- 
ilar ratio between cage size and particle size has been 
observed in (36j . However there is an important quali- 
tative difference between the motion of a particle in col- 
loidal suspension and other glass forming systems, and 
that observed in our system. In glass forming systems 
a particle spend most of its time rattling inside a cage. 
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FIG. 11: (Color online) Probability distribution of cos(a), 
where a is the angle formed by the displacements A of a 
particle in two subsequent taps. The distribution evidences 
that particles usually travel along straight lines, and that this 
tendency increases as the driving intensity decreases. 



from which it escape via infrequent cage-breaking rear- 
rangements. Here we observe the opposite behavior: par- 
ticles usually diffuse, and sometimes they get trapped in 
a cage. This unusual property of the trajectory is put in 
evidence by the study of the distribution of the angle a 
formed by the displacement A„ of a particle during tap 
n, and the displacement A„-|_i of the same particle dur- 
ing the following tap. Figure. [Tl] shows the distribution 
of 
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FIG. 12: (Color online) Probability distribution of vertical 
{Az, upper panel) and horizontal {Ah, lower panel) displace- 
ments of a grain during the 1st, 10th, 100th and 300th tap in 
semi-logarithmic (main panels) and linear (insets) scale. 



which is strongly peaked near 1. This is a clear indica- 
tion that particles usually move along straight lines (as 
also confirmed by Fig. llOp . an that cage motion is neg- 
ligible: for a particle rattling in a cage two consecutive 
displacements have opposite direction, and cos (a) ~ — 1. 

The absence of cage motion is due to the driving mech- 
anism. During a tap the system expands, cages break, 
and it is easier for the grains to move one with respect 
to the other. 



B. Particle displacements 

Here we investigate the evolution of the distribution 
of particle displacements during a single tap. After the 
application of a flow pulse to our system (and the follow- 
ing relaxation) a particle i, initially located in fi, will be 
in a new position position -I- A^, where Aj denotes its 
displacement. We examine below the probability that a 
particle makes a given displacement A. Due to presence 
of gravity, which breaks the up-down symmetry of the 
system, it is convenient to separate the vertical compo- 
nent of the displacement, Azi, from the horizontal ones. 



Axi and Ai/i. As Axi and Ayi have the same, even distri- 
bution, we have studied the evolution of the distribution 
of Ah = |Aa;| (|A?/|). In order to follow the dynamics 
of the system, in Fig. [T^ we have plotted the probabil- 
ity distribution of Az, Pz(Az), (upper panel) and of Ah, 
Ph{Ah), (lower panel) during tap number 1, 10, 100 and 
300. 

During the compaction process displacements with 
Az < are more probable than those with Az > 0. 
Therefore, as shown in Fig. [T^ (upper panel), Pz{Az) 
is asymmetrical. The asymmetry of the distribution de- 
creases as the system compactifies and after 300 taps, 
when stationarity is almost attained, Pz{Az) appears to 
be nearly symmetric. Accordingly the value of Az where 
Pz{Az) has its maximum increases, starting from a neg- 
ative value, until it reaches zero. Also, it is apparent 
that as the system compactifies the variance of the dis- 
tribution decreases. The probability of a large vertical 
displacement to occur is smaller in a dense rather than 
in a fluffy system. This is also true for the probability of 
large horizontal displacements, as shown in Fig. 1121 (lower 
panel) . 
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An important feature of the probability of both verti- 
cal and horizontal displacements is the nearly exponen- 
tial decay at large displacements. This is an indication 
of the fact that, while during a tap most of the parti- 
cles are subject to small displacements, very few of them 
may move much more. We conclude that the system is 
characterized by an heterogeneous dynamics. This is a 
well know property of thermal systems, like liquids or 
colloids (see [37| for a review) , which appears upon cool- 
ing the system near the glass transition. For instance, in 
colloidal systems the exponential tail of the particle dis- 
placement distribution has been experimentally observed 
by Weeks et. al. [H]. 

The probability distribution function of particles dis- 
placements can be further analyzed via the study of the 
excess- kurtosis 



' 3((A,-(A,))4) 



1, 



/i,z. 



(6) 



a comparison between the second and the fourth central 
moment of the distribution, which is zero for a gaussain 
distribution. We have first considered the excess-kurtosis 
as obtained from the probability of particle displacements 
in a single tap: a!^ and fluctuate from tap to tap, and 
their mean values are = 2.8±0.1, = 4.3±0.2. This 
positive values indicate that the probability distribution 
of particle displacements is more peaked with respect to 
a gaussian distribution. Then, we have considered the 
excess-kurtosis of the probability distribution of the hor- 
izontal and vertical components of 'ri(n)— ri(no), the total 
displacement of a grain after tap np, where uq — 100 cor- 
respond to the compaction time. As expected form the 
central limit theorem at long times this excess-kurtosis is 
zero: Fig. [T3l shows, interestingly, that in our system the 
excess kurtosis approach zero with a monotonic decay. 
This is in sharp contrast with the observations made in 
glass forming system, as in Ref . 36^ , where a peak is ob- 
served at the a relaxation time, i.e., when cage rearrange- 
ments occur. Therfore Fig. [T3] confirms the marginal role 
played by cage motion in our system. 

Our results put in evidence a very smooth behavior 
of the particles displacement distribution. Particularly 
we have not observed any 'rare event' (displacement of 
the order of 0.4 particle diameters), recently observed 
by Ribiere and coworkers [ssj in the study of a granular 
system undergoing compaction. This is probably due to 
the different driving mechanism of the systems. In shaken 
system, in fact, grains move one with respect to the other 
mainly when the pack settle downs and a shock wave 
propagates upwards; in our system, on the contrary, there 
is not a shock wave propagating as grains settle down 
slowly, and relative grain motion occurs during the tap. 



C. Spatial heterogeneous dynamics 

In supercooled liquids and dense colloidal systems the 
dynamics is heterogeneous as there are both slow and 




n-n„ 



FIG. 13: Time variation of the excess kurtosis of the probabil- 
ity distributions of the total horizontal (02 ) and vertical (of) 
displacements of a grain since tap no — 100. The monotonic 
decay evidences the marginal role played by cage motion in 
our system. 



fast particles. Moreover, fast particles are know to be 
spatially correlated as they appear to form clusters [s^, 
[33, |3^. Here we show that the exact same tendency 
characterizes a granular system subject to vertical taps, 
as suggested in Ref.s [s^ 0| . 

In order to characterize the spatially heterogeneous dy- 
namics [H, |4l|, li^l one usually resorts to a four-point 
time-dependent density correlation function and to its 
fiuctuations (susceptibility). This latter measures the 
correlated motion between pairs of particles. As this 
motion is decorrelated on short and long times, the sus- 
ceptibility show a well defined maximum at a given time. 
Unfortunately we cannot follow this line here as many av- 
erages are needed in order to get clear data on the fluctu- 
ations, and our simulations are computationally too ex- 
pensive. Moreover during compaction the system is not 
in a stationary state. In order to measure the degree of 
spatally heterogeneous dynamics we have therefore used 
a different method, based on the comparison between our 
system and a random one, as discussed below. 

We apply k flow pulses {k = 1,10,100,300) 
to our system (who reaches volume fraction 
0.576,0.596,0.612,0.616) and we determine the Up 
particle having experienced the largest displacement 
during the last tap. Here Up = pN, were N = 1600 is 
the total number of particles, p = 1%, 2%, 5% and 10%. 
When these fast particles are drawn, as in Fig. [141 it 
is apparent that they form clusters, an clear evidence 
of the spatially heterogeneous dynamics. We have 
quantified the degree of heterogeneity of the system as 
follows. After every tap we have determined the Up faster 
particles of the system, and determined the number 
Sp of couples of these particles made of neighboring 
particles (we consider two particles to be neighbor if 
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1% 2% 5% 10% 

FIG. 14: Granular compaction is characterized by a spatially heterogeneous dynamics. After the fcth tap (fc = 1, 10, 100,300, 
and volume fraction <1> = 0.576, 0.596, 0.612, 0.616) we plot the position of the rip particles which have experienced the largest 
displacement during the last tap, with rip = pN and p = 1%, 2%, 5% and 10%. 

the distance between their center is smaher than 1.2 quantity g^'^idom ^.j^g ^^^^ randomly selected 

particle diameters). Then we have computed the same particles of the system. A measure of the degree of 
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FIG. 15: (Color online) Granular compaction is characterized 
by a spatially heterogeneous dynamics. This is quantified by 
the parameter (see EqO which is plotted here as a function 
of the number of taps for p — 1%, 2%, 5% and 10% (from top 
to bottom). 



FIG. 16: (Color online) The mean displacement A of 
particle during a single tap increases with its Voronoi' 
volume V. The displacements are measured during tap 
1, 10, 100 and 300, when the volume fraction of the system 
is 0.576,0.596,0.612,0.616. 



spatial heterogeneous dynamics is given by 

„ ^random 

£ ^ fP (7) 

^P ^ ' 

Clearly ~ for an homogeneous systems, <^p > if 
faster particles form clusters, while < if faster par- 
ticles tend to be apart. Fig. [T5l shows the evolution of ^p 
with the number of taps. In all cases ^p > 0, signaling the 
presence of an heterogeneous dynamics. When the sys- 
tem approaches the steady state, and compaction stops, 
^p fluctuates around a plateau which varies with p be- 
tween 0.5 and 1.5. These are very high values, indicating 
a high degree of spatial heterogeneity of the system. 



D. Volume-Displacement correlation 

It is well known than many equilibrium and transport 
properties of dense fluids depend on the space available 
for molecular motion. For instance the well known free 
volume theory developed by Choen and TurnbuU [43] to 
explain the divergence (with a Vogel-Tamman-Fulcher 
law) of the relaxation time of many glass formers as the 
temperature is decreased, is based on the idea that the 
space available for molecular motion decreases with the 
temperature. It is therefore interesting to check for cor- 
relations between the displacement of a particle and its 
free volume in our granular system. There arc many 
possible way to define the free volume of a particle. Here 
we approximate, for simplicity sake, we consider the free 
volume of particle i to be V/" = Vi — Vq, where Vi is 
the Voronoi volume of the particle (see Sec. lIVBp be- 
fore the application of a tap, and Vq = 4/37r(£'/2)'^ 
its volume. The displacement A; of particle i (A^ = 



(A^ -I- A^ -I- A^)^/^) is the distance between the position 
of the particle before and after the application of the tap. 

A possible connection between Vf' and A^ is suggested 
by the similarity between the the probability distribution 
functions of Voronoi volumes (sec Fig. [7la), and that of 
particles displacemets (Fig. [T^ . Both of them have an 
exponential tail at high values. Moreover they evolve in 
a qualitatively similar way (the variance and the mean 
value decrease) as the system compactifies. 

In order to test this possible correlation we have com- 
puted during the nth tap {n = 1, 10, 100, 300) (tq = 0.03 
s, V = 0.2 cm/s) the mean value of the displacement 
A of all the particles with free volume . The depen- 
dence of A on is shown in Fig. 1161 This figure puts in 
evidence the existence of an almost linear correlation be- 
tween Voronoi volumes and particle displacements, with 
dA/dv ~ 0.01/0.2 = 0.05 cm-^: the larger the Voronoi' 
volume of a particle, the bigger its displacement. 



E. Force-Displacement correlation 

In the previous paragraph we have shown that there 
is a correlation between the displacement A of a particle 
during a single tap, and its free volume. Fig. [T7] inves- 
tigates the correlation between displacement and com- 
pression C (see Eq. S]) of a particle. The figure puts in 
evidence that for small values of the compressional force 
there is a decreasing linear relation between displacement 
and compression of a particle, while for higher value of 
the compressional force the two becomes uncorrelated. 
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FIG. 17: (Color online) The mean displacement A of a grain 
as a function of its compressional force. When the packing is 
loose grains with higher compressional forces move less. As 
the granular pack compactifies the displacement of a particle 
appears not to depend on its compression. 



VI. CONCLUSIONS 

We have reported results of a numerical simulations 
of a packing of monosize spheres submitted to vertical 
taps made of flow pulses as in the experiment of Ref. @ . 
Our results relative to the dynamics of the systems con- 
firms earlier experimental observation: as the intensity 
of vibration decreases both the value of the volume frac- 
tion reached at stationary and the compaction time in- 
creases [l], H, H, d, [3| • The increase of the compaction 
time with the decreasing of the vibration intensity is dra- 
matic as this appears to diverge with a power law when 
the fluid velocity goes to zero. 

The analysis of the evolution of several structural 
quantities during compaction has revealed that this is not 
accompanied by any particular geometrical modification. 
The radial distribution function and the Voronoi volume 
distribution, for instance, smoothly change as the den- 



sity of the granular system increases. In particular the 
collapse of the Voronoi volume distributions (sec. lIVBp 
evidences the presence of a single underlying geometri- 
cal structure in the system. Also, the probability distri- 
bution function of normal forces between grains do not 
change during compaction. 

The analysis of the dynamics of compaction has re- 
vealed that this is characterized by dynamical hetero- 
geneities. The probability that, during a tap, a particle 
makes a given displacement decreases exponentially with 
its size, resembling observations made in dense colloidal 
systems [36| . Similarly we have observed that, during a 
tap, faster particle tend to form cluster, as observed both 
in colloidal [s^ and in glass-forming [33,[3l,|4l| systems. 
There is, however, a marked difference between a typical 
trajectory of a granular particle during compaction, and 
a typical trajectory of a particle in supercooled liquids. 
Precisely, in supercooled liquids a particle spend most 
of time in cages formed by its neighbors, and occasion- 
ally makes large displacement escaping from the cage. In 
our system, instead, particles usually diffuse, and some- 
times are trapped in a cage. This different behaviour is 
due to the particular driving of our system as in our sys- 
tem, when the the flow in on the system expands and the 
grains are able to escape from their cages. 

As the slowdown of the dynamics is related to the space 
available for particle motion, we have studied the corre- 
lation between the displacement of particle during a tap 
and its Voronoi volume, which is a rough estimate of its 
free volume. This analysis has shown that particles with 
larger Voronoi' volumes are those who make larger dis- 
placements during a tap. 
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